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Abstract 

The existence of dark matter (DM) at scales of few pc down to ~ 10^^ pc around the centers 
of galaxies and in particular in the Galactic Center region has been considered in the literature. 
Under the assumption that such a DM clump, principally constituted by non-baryonic matter 
(like WIMPs) does exist at the center of our galaxy, the study of the 7-ray emission from the 
Galactic Center region allows us to constrain both the mass and the size of this DM sphere. 
Further constraints on the DM distribution parameters may be derived by observations of bright 



used estimates of the enclosed 



infrared stars around the Galactic Center. Hall and Gondolo 
mass obtained in various ways and tabulated by Chez et al. l3|. Moreover, if a DM cusp does 
exist around the Galactic Center it could modify the trajectories of stars moving around it in a 
sensible way depending on the DM mass distribution. Here, we discuss the constraints that can 
be obtained with the orbit analysis of stars (as S2 and S16) moving inside the DM concentration 
with present and next generations of large telescopes. In particular, consideration of the S2 star 
apoastron shift may allow improving limits on the DM mass and size. 
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I. INTRODUCTION 



In the last years intensive searches for dark matter (DM), especially its non-baryonic 
component, both in galactic halos and at galaxy centers have been undertaken (see for 
example {4, 5 1 for recent results). It is generally accepted that the most promising candidate 
for the DM non-baryonic component is neutralino. In this case, the 7-flux from galactic halos 
and from our Galactic halo in particular) could be explained by neutralino annihilation 
, lid . I 111 . I12I. Il3[|. Since 7-rays are detected not only from high gala ctic latitude, but 



also from the Galactic Center, there is a wide spread hypothesis (see IJ] for a discussion) 
that a DM concentration might be present at the Galactic Center. In this case the Galactic 



Center could be a strong source of 7-rays and neutrinos 



fl 0, H, Q, 0, Q, Q: 
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due to DM annihilation. Since it is also expected that DM forms spikes at galaxy centers 



23|, 



24 



25l | the 7-ray flux from the Galactic Center should increase significantly in that case. 
At the same time, prog ress in monitoring bright stars near the Galactic Center have been 



reached recently 



ime, progi 



The astrometric limit for bright stellar sources near the Galactic 



Center with 10 meter telescopes is today 66iq ~ 1 mas and the Next Generation Large 



Te 



escope (NGLT) will be able to improve this number at least down to ^^30 ~ 0.5 mas 



28l . |29[ | or even to ^^30 ~ 0.1 mas 27], |28|, |29(] in the K-band. Therefore, it will be possible 



to measure the proper motion for about ~ 100 stars with astrometric errors several times 
smaller than errors in current observations. 

The aim of this paper is to constrain the parameters of the DM distribution possible 
present around the Galactic Center by considering the induced apoastron shift due to the 
presence of this DM sphere and either available data obtained with the present generation 
of telescopes (the so called conservative limit) and also expectations from future NGLT 
observations or with other advanced observational facilities. 



II. THE MASS CONCENTRATION AT THE GALACTIC CENTER 

Recent advancements in infrared astronomy are allowing to test the scale of the mass 
profile at the center of our galaxy down to tens of AU. With the Keck 10 m telescope, the 
proper motion of several stars orbiting the Galactic Center black hole have been monitored 
and almost entire orbits, as for example that of the S2 star, have been measured allowing an 
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unprecedent description of the Galactic Center region. Measurements of the amount of mass 
M(< r) contained within a distance r from the Galactic Center are continuously improved 
as more precise data are collected. Recent observations [2] extend down to the periastron 
distance (~ 3 x 10~^ pc) of the S16 star and they correspond to a value of the enclosed mass 
within ~ 3 X 10^'* pc of ^ 3.67 x 10^ Mq. Several authors have used these observations to 
model the Galactic Center mass concentration. Here and in the following, we use the three 
component model for the central region of ourgalaxy based on estimates of enclosed mass 
given by Chez et al recently proposed This model is constituted by the central 

black hole, the central stellar cluster and the DM sphere (made of WIMPs), i.e. 



M(< r) = Mbh + M,(< r) + MnMi< r) , 



(1) 



where Mbh is the mass of the central black hole Sagittarius A*. For the central stellar 
cluster, the empirical mass profile is 



MJ< r) 



1.6 



1.0 



r > 



(2) 



with a total stellar mass = 0.88 x 10^ Mq and a size R^: = 0.3878 pc. 

As far as the mass profile of the DM concentration is concerned. Hall and Gondolo |l| 
have assumed a mass distribution of the form 



Mdm{< r) 



DM 



Rdm 



3-a 



r <R 



DM 



(3) 



DM, 



r > R 



DM 



Mum and Rdm being the total amount of DM in the form of WIMPs and the radius of the 
spherical mass distribution, respectively. 

Hall and Gondolo l| discussed limits on DM mass around the black hole at the Galactic 
Center. It is clear that present observations of stars around the Galactic Center do not 
exclude the existence of a DM sphere with mass ~ 4 x 10^ Mq, well contained within the 
orbits of the known stars, if its radius Rdm is ^ 2 x 10""^ pc (the periastron distance of 
the S16 star in the more recent analysis |3|). However, if one considers a DM sphere with 
larger radius, the corresponding upper value for M^m decreases (although it tends again to 
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increase for extremely extended DM configurations with Rdm ^ 10 pc). In the following, 
we will assume for definiteness a DM mass Mom ~ 2 x lO^M©, that is the upper value for 
the DM sphere in [l] within an acceptable confidence level in the range 10^'^ — 10^^ pc for 
Rdm- As it will be clear in the following, we emphasize that even a such small value for the 
DM mass (that is about only 5% of the standard estimate 3.67 ±0.19 x 10^ Mq for the dark 
mass at the Galactic Center jsj) may give some observational signatures. 



Evaluating the S2 apoastron shift 



38| function of Rdm, one can further constrain 



the DM sphere radius since even now we can say that there is no evidence for negative 
apoastron shift for the S2 star orbit at the level of about 10 mas 26|]. In addition, since at 
present the precision of the S2 orbit reconstruction is about 1 mas, we can say that even 
without future upgrades of the observational facilities and simply monitoring the S2 orbit, 
it will be possible within about 15 years to get much more severe constraints on Rdm- 
Moreover, observational facilities will allow in the next future to monitor faint infrared 



objects at the astrometric precision of about 10 /ias |30| and, in this case, previous estimates 
will be sensibly improved since it is naturally expected to monitor eccentric orbits for faint 
infrared stars closer to the Galactic Center with respect to the S2 star. 

In Fig. [H the mass profile M(< r) obtained by using observations of stars nearby the 
Galactic Center is shown (solid line). The dotted line represents the stellar mass profile as 
given in Eq. ([2]), while the dashed lines are for DM spheres with mass Mdm — 2 x 10^ Mq 
and radii Rdm = 10^'^ and 10^^ pc, respectively. 

In the following section, we study the motion of stars as a consequence of the gravitational 
potential $(r) due the mass profile given in Eq. ([T]). As usual, the gravitational potential 
can be evaluated as 

^{r) = -G j^^dr' . (4) 

For convenience, in Fig. [2] the gravitational potential due to the total mass (solid line) 
contained within r is given as function of the galactocentric distance. For comparison, the 
contributions due to the single mass components, i.e. central black hole, stellar cluster and 
DM sphere, are also shown. 
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FIG. 1: The mass M(< r) obtained in [2'] from observations of stars at the Galactic Center is 
shown (sohd hne). The dotted hne represents the stellar mass profile as given in Eq. ([2]), while the 
dashed lines are for DM spheres with radii Rdm = 10"'^ and 10~^ pc and mass Mjjm — 2 x 10^ 
Mq, that corresponds to some acceptable estimate for the upper limit of Mom from Figs. 4-6 in 
l|] for Rdm in the above range of values. 

III. APOASTRON SHIFT CONSTRAINTS 



According to GR, the motion of a test particle can be fully described by solving the 
geodesic equations. Under the assumption that the matter distribution is static and pres- 
sureless, the equations of motion at the first post-Newtonian (PN) approximation become 



[see e.g. 
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36|) 



^ ~ -V($7v + 2$^) + 4v(v ■ V)$7v - v^V^N 



(5) 



We note that the PN-approximation is the first relativistic correction from which the apoas- 



tron advance phenomenon arises 



39|. In the case of the S2 star, the apoastron shift as 



seen from Earth (from Eq. ([7])) due to the presence of a central black hole is about 1 mas, 
therefore not directly detectable at present since the available precision in the apoastron 
shift is about 10 mas (but it will become about 1 mas in 10-15 years even without consid- 
ering possible technological improvements). It is also evident that higher order relativistic 
corrections to the S2 apoastron shift are even smaller and therefore may be neglected at 
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FIG. 2: The gravitational potential $(r) (solid line) in cgs units as a function of the galactocentric 
distance r as due to the mass M{r) in Eq.([T]) is shown. For comparison, also the gravitational 
potentials due to the single mass components, i.e. black hole (dashed line), stellar cluster (dot- 
dashed line) and DM (dotted line), are also given. Here we assume that DM mass Mdm — 2 x 10^ 
M© and radius Rdm = pc. 



present, although they may become important in the future. 

As it will be discussed below, the Newtonian effect due to the existence of a sufficiently 
extended DM sphere around the black hole may cause a apoastron shift in the opposite 
direction with respect to the relativistic advance due to the black hole. Therefore, we have 
considered the two effects comparing only the leading terms. ^ 

For the DM distribution at the Galactic Center we follow Eq. ([3]) as done by [L]. Clearly, 
if in the future faint infrared stars (or spots) closer to the black hole with respect to the 
S2 star will be monitored SOj, this simplified model might well not hold and higher order 
relativistic corrections may become necessary. 

For a spherically symmetric mass distribution (such as that described above) and for a 
gravitational potential given by Eq. (jlj), Eq. ([5]) may be rewritten in the form (see for details 

m 



dv 



GM{r) 



1 H — H 



4v(v • r) 



(6) 



r and v being the vector radius of the test particle with respect to the center of the stellar 
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cluster and the velocity vector, respectively. Once the initial conditions for the star distance 
and velocity are given, the rosetta shaped orbit followed by a test particle can be found by 
numerically solving the set of ordinary differential equations in eq. (jS])- 

In Fig. [31 as an example, assuming that the test particle orbiting the Galactic Center 
region is the S2 star, we show the Post Newtonian orbits obtained by the black hole only, 
the black hole plus the stellar cluster and the contribution of two different DM mass density 
profiles. In each case the S2 orbit apoastron shift is given. As one can see, for selected 
parameters for DM and stellar cluster masses and radii the effect of the stellar cluster 
is almost negligible while the effect of the DM distribution is crucial since it enormously 
overcome the shift due to the relativistic precession. Moreover, as expected, its contribution 
is opposite in sign with respect to that of the black hole 33|. 

We note that the expected apoastron (or, equivalently, periastron) shifts 
(mas/revolution), A$ (as seen from the center) and the corresponding values as seen 
from Earth (at the distance Ro — 8 kpc from the GC) are related by 

A0| = i^lMA*. (7) 

Kq 

where with the sign ± are indicated the shift angles of the apoastron (+) and periastron 
respectively. The S2 star semi-major axis and eccentricity are d = 919 AU and e = 0.87 

In Fig. m the S2 apoastron shift as a function of the DM distribution size Rdm is given 
for a = and M^m — 2 x 10^ Mq. Taking into account that the present day precision for 
the apoastron shift measurements is of about 10 mas, one can say that the S2 apoastron 
shift cannot be larger than 10 mas. Therefore, any DM configuration that gives a total 
S2 apoastron shift larger than 10 mas (in the opposite direction due to the DM sphere) is 
excluded. The same analysis is shown in Figs. [5] and [6] for two different values of the DM 
mass distribution slope, i.e. a = 1 and a = 2, respectively. In any case, we have calculated 
the apoastron shift for the S2 star orbit assuming a total DM mass M^m — 2 x 10^ Mq. As 
one can see by inspecting Figs. HJIHI the upper limit of about 10 mas on the S2 apoastron 
shift may allow to conclude that DM radii in the range about 10^^ — 10^^ pc are excluded 
by present observations. 

We notice that the results of the present analysis allows to further constrain the results 
of the Hall and Gondolo paper [l| who have concluded that if the DM sphere radius is in the 
range 10^^ — 1 pc, configurations with DM mass up to Mdm = 2 x 10^ Mq are acceptable. 
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The present analysis shows that DM configurations of the same mass are acceptable only 
for Rdm out the range between 10~^ — 10~^ pc, almost irrespectively of the a value. 



IV. CONCLUSIONS 



In this paper we have considered the constraints that the upper limit (presently of about 
10 mas) of the S2 apoastron shift may put on the DM configurations at the galactic center 
considered by Hall and Gondolo 

When (in about 10-15 years, even without considering improvements in observational 
facilities) the precision of S2 apoastron shift will be about 1 mas (that is equal to the 
present accuracy in the S2 orbit reconstruction) our analysis will allow to further constrain 
the DM distribution parameters. In particular, the asymmetric shape of the curves in Figs. 
llUni imply that any improvement in the apoastron shift measurements will allow to extend 
the forbidden region especially for the upper limit for Rdm- 

In this context, future facilities for astro metric measurements at a level 10 /ias of faint 
infrared stars will be extremely useful jsO] and they give a chance to put even more severe 
constraints on DM distribution. 

In addition, it is also expected to detect faint infrared stars or even hot spots 3J] orbiting 
the Galactic Center. In this case, consideration of higher order relativistic corrections for 
an adequate analysis of the stellar orbital motion have to be taken into account. 

In our considerations we adopted simple analytical expression and rehable values for 
Rdm and Mdm parameters following [l| just to illustrate the relevance of the apoastron 
shift phenomenon in constraining the DM mass distribution at the Galactic Center. If other 



models for the DM distributions are considered (see, for instance [35| and references therein) 
the qualitative aspects of the problem are preserved although, of course, quantitative results 
on apoastron shifts may be different. 
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FIG. 3: PN-orbits for different mass configurations at the Galactic Center. The S2 star has been 
considered as a test particle and its apoastron shift is indicated in each panel as A<1> (in arcsec). 
The top-left panel shows the central black hole contribution to the S2 shift that amounts to about 
580 arcsec. The top-right panels shows the combined contribution of the black hole and the stellar 
cluster (taken following eq. [2]) to the S2 apoastron shift. In the two bottom panels the contribution 
due to two different DM mass-density profiles is added (as derived in eq. [3]). We assume that DM 
mass Mdm ^ 2 X 10^ Mq. 



10 




0.001 0.010 0.100 1.000 10.000 

Rdm (pc) 



FIG. 4: Apoastron shift as a function of the DM radius Rdm for a = and M^m — 2 x 10^ Mq. 
Taking into account present day precision for the apoastron shift measurements (about 10 mas) 
one can say that DM radii Rdm in the range 8 x 10^'^ — 10~^ pc are not acceptable. 




FIG. 5: The same as in Fig. [J] for a = 1 and Mdm — 2 x 10^ Mq. As in the previous case one can 
say that the S2 apoastron shift put severe limits on the DM mass radii that are not acceptable in 
the range 9 x 10~^ — 10^^ pc. 
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FIG. 6: The same as in Fig. [3]for a = 2 and Mj^m ~ 2 x 10^ Mq. As in the previous case one can 
say that the upper hmit to the S2 apoastron shift allows to constrain the DM radius to be out the 
range 1.0 x 10^'^ — 1.1 x 10^^ pc. 
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